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Abstract: It is a widely observed phenomenon in computer graphics that the size of the 
silhouette of a polyhedron is much smaller than the size of the whole polyhedron. This paper 
provides, for the first time, theoretical evidence supporting this for a large class of objects, 
namely for polyhedra that approximate surfaces in some reasonable way; the surfaces may 
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polyhedra have silhouettes of expected size O(y^) where the average is taken over all points 
of view and n is the complexity of the polyhedron. 
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Une borne suprieure sur la taille moyenne des 

silhouettes 

Resume : II est connu en infographie que la taille de la silhouette d'un polydre s'avre 

soiivcnt, cn pratique, bicn plus petite que la taille du polydre entier. Get article est le 
premier fournir des preuves thoriques justifiant cette observation pour une large classe 
d'objets : les polydres qui approximent des surfaces de manire raisonnable ; les surfaces 
considres ne sont pas nccssairement convexcs ou lisscs ct dies pcuvent avoir un bord. Nous 
prouvons qu'un tel polydre de taille n a, en moyennant sur tous les points de vue I'infini, 
une silhouette de taille 0{^/n). 

Mots-cles : silhouette, polydre, borne suprieure 



An Upper Bound on the Average Size of Silhouettes 



3 



1 Introduction 

The silhouette of a polyhedron with respect to a given viewpoint is, roughly speaking, the 
set of edges incident to a front and a back face. Silhouettes arise in various problems in 
computer graphics such as hidden surface removal and shadow computations (see [HI HI [S] for 
some recent references) and algorithms to compute them efficiently have been well-studied 
(see the survey by Isenberg et al. [6]). They are important in shape recognition; Sander 
et al. [9] claim that the silhouette "is one of the strongest visual cues of the shape of an 
object" . 

It is a widely accepted fact that the silhouette of a polyhedron is usually much smaller 
than the whole polyhedron. Sander et al. [3], for instance, state the largely repeated claim 
that the silhouette of a mesh is often of size ld{y/n) where n is the number of faces of the 
mesh. An experimental study by Kettner and Welzl [7] confirms this for a set of realistic 
objects. This experimental study was extended by McGuire [8] to a larger database of larger 
objects for which the observed size of the silhouette is approximately 

There are few theoretical results supporting these observations. Kettner and Welzl [7] 
prove that a convex polyhedron that approximates a sphere with Hausdorff distance e has 
<d{l/e) edges, and a random orthographic projection of such a polytope has Q{l/^/e) sil- 
houette edges. Alt et al. [1] give conditions under which it can be proved that the average 
silhouette of a convex polyhedron has size 0{^/n) and give additional conditions under which 
the worst-case size is provably sub-linear. 

The goal of this paper is to study the average silhouette size of non- convex polyhedra. 
Convexity is a very strong assumption, which was crucial in the previous theoretical results. 
Here, rather, we assume that the polyhedron is a good approximation of some fixed (not 
necessarily convex) surface. Notice that it is very difficult to guarantee anything on the 
worst- case complexity of the silhouette of a polyhedron unless it approximates a strictly 
convex surface. Alt et al. \T\ give an example of a polyhedral approximation of a section of 
a cylinder with worst-case silhouette size Q{n) (see Figure [1]). Moreover, their example can 



Figure 1: A worst-case linear silhouette (left) of a polyhedron approximating a cylinder. 
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be modified in such a way that the surface is smooth, and its polyhedral approximation is 
as "nice" as one might hope (for instance, it can be required that the faces are fat triangles 
that all have almost the same size). 

In this paper we prove an upper bound on the expected size of the silhouette for random 
viewpoints. We prove that the silhouette of a polyhedron that approximates a surface in a 
reasonable way has expected size 0{-\/n). Note that the average is taken over all viewpoints 
for a given surface, and not on a set of surfaces. 

In Section [21 we define precisely the notion of silhouette for polyhedra and general sur- 
faces. We then present and prove our main result in Section [3] and conclude in Section |4l 



2 Definitions 

The term silhouette has been used in the literature to represent several different notions, 
depending on the application, refiecting such issues as: is the object considered opaque or 
transparent? Is occlusion taken into account? Is one interested by what the eye perceives, 
i.e., a plane curve, or by the space curve which gave birth to it? In the area of photography, 
for instance, a silhouette (also called apparent boundary) is defined as an outline of a solid 



object, as cast by its shadow, that appears dark against a fight background (Figure 2(a) I. 
In the field of computer vision, by contrast, the silhouette (also called rim, profile or contour 
generator) is roughly defined as the curve on the surface that separates front face regions 



from the back ones, either for opaque (Figure 2(b) I or for transparent (Figure 2(c) I objects 



In this paper we prove an upper bound on the size of the transparent silhouette; since 
such a silhouette contains the apparent boundary and the contour, our bounds also apply 
to all these types of silhouettes. In the rest of the paper the term silhouette will be used to 
mean transparent silhouette. 
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In the rest of this section we give a formal definition of silhouettes of polyhedra and then 
provide a definition for more general surfaces. 

2.1 Polyhedra 

The (transparent) silhouette of a polyhedron from a viewpoint (possibly at infinity) is the 
set of edges that arc adjacent to a front face and a back face. A face is considered a front 
face if the angle between its normal vector and a vector from a point of the face to the 
viewpoint is acute, and a back face if that angle is larger than 7r/2. If the point of view is in 
the plane containing the face, we refer to the definition of silhouettes for the case of general 
surfaces. The normal vectors should point outwards, but what really matters is that the 
orientation is consistent for the two faces that share this edge, so this definition also applies 
to non-orientable (necessarily self-intersecting) polyhedra. 

In this paper, we call complexity of a silhouette (of a polyhedron) its number of edges. 

2.2 General surfaces 

Our objective is to bound the size of the silhouette of a polyhedron. To achieve this goal, 
we need to relate the silhouette of the polyhedron to the silhouette of the surface it ap- 
proximates, which means we need a definition of silhouettes that applies to a larger class of 
objects. Although this may seem unintuitive, we first define the silhouette as a set of rays, 
and then relate this to the more usual concept of a set of points on the surface. 

Let 5 be a compact 2-manifold without boundary. It separates M.^ in two non-empty 
open regions; call and ^' their closures (so ^ n ^' = 5 and ^ U ff' = K^). Let V be 
a viewpoint not on S but possibly at infinity. The (transparent) silhouette of S from V 
is the set of rays R starting from V that are tangent to 5* in a non-crossing way {R may 
cross S elsewhere). More formally, we require that there exists an open segment u of R that 
contains a connected component of i? n 5 and is contained either in ff or . 

This definition defines a set of rays. The silhouette can also be seen as the trace of this 
set of rays on the surface. More precisely, for each ray R on the silhouette, we consider 
the closest point to V on each connected component oi RO S that satisfies the non-crossing 
property. This definition is consistent with the one given for the particular case of polyhedra, 
and is the one we will use in this paper. 

For a given viewpoint at infinity, we define the (projected) length of the silhouette as the 
length (counted with multiplicity if several points have the same projection) of the projection 
of the silhouette, along the direction given by the viewpoint, on an orthogonal plane. 

Remark. The definition of the silhouette can be extended to cases where S is not a 2- 
manifold, but an immersion of a compact 2-manifold. More precisely, we have a 2-manifold 
S' and an application f : S' ^ M.^ such that S = f{S') and for any point on S' there 
exists a neighborhood U of that point such that U and f{U) are homeomorphic. The local 
orientation is sufficient to decide whether R crosses S or not (note that more complicated 
things can happen than crossing or being tangent, even with smooth surfaces; for instance. 
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the surface may ripple an infinite number of times in the neighborhood of a point, making 
it impossible to define on which side of S' i? is near the intersection point). This remark 
extends to the whole paper and, in particular, to Theorem[T] However, we do not give either 
a definition or a proof of this, as it would uselessly make everything more obscure. 

3 Main results 

Let S* be a compact 2-manifold without boundary whose silhouettes have finite average 
length, silh(5), where the average is taken over all viewpoints at infinity. Let P„ be a 
polyhedron with n triangular faces, that is homeomorphic to S through /„ : P„ S, such 
that: 

1. the length of any edge of P„ is at least and 

2. for any point a; on P„, d{x,fn{x)) < ^-^^ where h{x) is the smallest height of the 
triangle(s) of Pn that contain(s) x, 

where a and /? are two arbitrary positive numbers and d() denotes the Euclidean distance. 

Theorem 1. The expected complexity of the silhouette of Pn is 0{y/n), where the average 
is taken over all viewpoints at infinity. More precisely, for any n, the expected complexity is 
at most 



Note that the bound is valid for any n and any polyhedron P„ satisfying the above 
assumptions. Note also that the bound depends on S only by the average length of its 
silhouette. 

We first clarify the meaning of the hypotheses on P„ and their implications. We then 
prove Theorem [1] in Section 13.21 We finally show in Section 13.31 how Theorem [1] can be 
generalized to surfaces with boundary and viewpoints at finite distance. In particular, we 
prove the following result. 

Let S" be any compact two-manifold with boundary of finite length and whose silhouette 
has finite average length (taken over all viewpoints at infinity). 

Theorem 2. Any mesh Pn with n triangular faces that approximates S" according to Hy- 
potheses 1 and 2 has a silhouette of expected complexity 0{y/n) when the viewpoint is chosen 
uniformly at random in a hall. 

3.1 Meaning of the hypotheses 

Hypothesis[T]is here to avoid short edges. The main idea of the proof is to link the complexity 
of the silhouette to its length, and arbitrarily short edges would make this impossible. Now 
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the factor makes sense: intuitively, since the polyhedron has n faces, each face has area 
of order - , which means that the edges have length of order . 

Hypothesis [2] is rather technical, and we discuss instead the meaning of the following two 
more intuitive hypotheses, which, together with Hypothesis [TJ impljQ Hypothesis [3 

3. The faces of P„ are fat. 

4. For any x on P„, d{x, fn{x)) < ^, where 7 is some positive constant. 

Hypothesis [3] is quite natural. Hypothesis [4] ensures that P„ approximates S. Further- 
more, the factor is reasonable; indeed, in 2D, when considering a regular polygon with 
edge length 0(-^) inscribed in a circle of radius 1, the maximal distance between a point 

on the polygon and the circle is 6(^)- The situation is the same in 3D. Basically it means 
that the error when approximating the surface with a plane is of the second order. 

Our hypotheses ([l][3]|4] or [TJll) ensure that the homeomorphism /„ has good properties, 
that is that, roughly speaking, the polyhedron can be obtained by only a small perturbation 
of the surface while keeping the normal vectors in approximately the same directions. This 
is crucial for our proof since otherwise, for example, a cylinder can be approximated by a 



lantern of Schwarz [TD] (see Figure 3(a) ) whose silhouette has expected complexity Q{n) and 
unbounded length. 

Notice that the existence of polyhedra with arbitrarily large number of edges that ap- 
proximate the surface according to these hypotheses is a constraint on the surface. Not 
every surface admits such an approximation (think of the neighborhood of in the surface 



^Indeed, for any x in P„, Hypotheses [T] and |3] imply that h{x) ^ S/\/n for some positive constant 5; 
Hypothesis |2] then follows from Hypothesis l4l since h{x)/^/n ^ S/n 5/7 ■ d{x, fn{x)). 
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Figure 4: (a) Length and dihedral angle of an edge; (b) set of directions for which e is on 
the silhouette. 



defined by z = {x^ + y^Y^^ as shown in Figure [3(b)[ ). However, the class of surfaces for 
which such approximations exist is quite large. It includes, in particular, smooth surfaces 
and polyhedra with fat faces. 



3.2 Proof of Theorem [T] 

We consider a point of view chosen randomly at infinity. We call le the length of an edge e 
of polyhedron P„ and 6e the exterior dihedral angle associated to e (see Figure |4]). 

Let Tg denote the union of the two triangles adjacent to edge e (including e but not 
the other edges). For any part 31 of 5*, let silh(IR) be the average length of the part of the 
silhouette of S that lies in 3?. 

We first recall a classical formula on the expected size of silhouettes which can also be 
found, for instance, in [8J. 

An edge e is on the silhouette if the direction of view is in the dark area of the sphere 
of directions of Figure [4(b)[ The angular measure of this region is AOe, which means that 
the probability for e to be on the silhouette is Oe/i^. The expected number of edges on the 
silhouette is thus 

TT 

edge e 

We now state our main lemma. The general idea of this lemma is that under strong 
hypotheses {S has bounded curvature, the edges have length 0(--i=), and Hypotheses [3] and 

m are satisfied), one can prove that 9^ ^ for some constant C. In cases where this 
inequality does not hold, edge e is near some kind of edge of the surface, or at least some 
feature that will appear quite often on the silhouette and we are going to charge this edge 
to the silhouette of S. 
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Figure 5: Construction of fl. 



Lemma 3. For any edge e on P„, 



^ ^ + — Sllh{fniTe)) with C= 31.3/3. 



Theorem [T] follows from Lemma|31 Indeed, since P„ has ^ edges, each of length at least 
(by Hypothesis 1), we get that the expected complexity of the silhouette is 

i;^i^^+8^3silh(5), 
TT 2 ^Jn a 

because ^ silh (/„ (Te)) ^ 3 silh(S') since the length of the silhouette of S that lies in the 

edge e 

image (through /„) of a triangle is counted three times (once per edge). Hence, 



Ei^ ( 15 /3 + — silh(5) ) ^=0 (V^) 



Proof of Lemma[3[ The idea of the proof is as follows. Consider the set of directions for 
which e is on the silhouette. We first construct a subset O of these directions whose measure 
is a constant times 9e — (see Figure [S]) . We then prove a lower bound on the length of 
the silhouette of fn{Te) for all these directions, and deduce the result. 

Let C be a positive constant, whose value will be defined later (see Equation [21) . For 
any edge e on P„, we can assume that 9^ — > since, otherwise, 6^ ^ and there is 
nothing else to prove. 
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The set of directions for which e is on the silhouette is the set of directions between the 
planes defined by the faces adjacent to e. Rotate each face about e by an angle of 2^ so 

that the exterior dihedral angle decreases by (see Figure 5(a) I. is defined to be the 



set of directions between these two new planes that make an angle larger than 7r/3 with 



the line supporting e; Figure 5(b) shows one component of O, the other one consists of the 
symmetric set of opposite directions. The measure of the set of directions between these two 
planes is A{9e — -^)- Restricting this set of directions to those that make an angle larger 
than 7r/3 with the line supporting e, we get, by integrating on the sphere of directions, that 
the measure of 57 is 2 (Oe — "^)- 

The remaining step uses the property, which we prove in Corollary [5l that for all the 
directions in fi, the silhouette of fn{Te) has length at least Ze/4. Assuming this temporarily, 
we sum this inequality over Vl. The smaller side of the inequality is 2^(^e — -^)- The 
larger side is the integral of the length of the silhouette of fn{Te) over all directions in Vi, 
which is smaller than this same integral over all directions, that is 47rsilh(/ji(Te)). Hence 
47rsilh(/„(Te)) ^ ^(^'e — -^), which concludes the proof. □ 

We now state a lemma and its corollary which we used in the proof of Lemma [3] under 
the hypothesis that 9e — > 0. We can thus assume in the sequel that this property holds. 

Let e' be the segment obtained by clipping from e all the points at distance less than ^■ 



from its extremities. Refer now to Figures 6(a) (b) 



Lemma 4. Any line with direction G that intersects e' can be translated in a direction 
orthogonal to e and d until it becomes tangent to S in fn{Te). 

Corollary 5. For any direction d in f2, the silhouette of fn{Te) has length at least ^. 
Proof. Consider the projection of e' and of the silhouette of fn{Te) onto a plane orthogonal 



to d (see Figure 6(a) I. It follows from Lemma |4] that, in that plane, each point on the 
projection of e' maps to a point on the projected silhouette in the direction orthogonal to 
e'. Hence, the projected silhouette is longer than the projection of e', which is at least 



times the length of e' since d makes an angle of at least 7r/3 with e'. Thus the silhouette of 

\/3k ^ k 
2 2^4 



/„(Te) has length at least > !f. □ 



Proof of Lemma^ Let D denote a line with direction d G fl that intersects e'. Let Ti 
and T2 denote the two triangles adjacent to e and let hi and /12 denote their respective 
smallest heights. Let Xi = Phi/y/n, = niax(xi,X2), and x^ — ™ii(xiiX2). Refer now 
to Figure [6(b)[ We call Dt, t G X^,X^], the line obtained by translating D at distance 
\t\ in a direction orthogonal to the plane defined by e and d; positive values of t correspond 
to lines in the half-space bounded by the plane defined by e and D, and not containing Tg, 
negative values of t correspond to lines in the other half-space. For clarity, we denote D-^^ 
by D- and by D+ . 

By construction, £)+ is at distance x~^ from T^. Thus £)+ does not intersect /„(Te), 
by Hypothesis 2. We prove that intersects fn{Te) and that no line Dt intersects the 
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Figure 6: For the proofs of Lemma H] and Corollary O 



boundary of fn{Te)- This will imply that, sweeping Dt from D+ to D~, the first line Dtg 
that intersects /„(Te) is tangent to /„(Te) at one of its interior point, which will conclude 
the proof. 

We first prove that no line Dt intersects the boundary of fn(Te). In other words, we 
prove that, for each edge e* on the boundary of Te, no line Dt intersects fn{e*)- Let Ti be 
the triangle (of Te) containing e*. By Hypothesis 2, it is sufficient to prove that the distance 
between Dt and e* remains greater than or equal to Xi for all t. 

First notice that it is sufficient to prove that the distance between Dt and e* remains 
greater than or equal to Xi for t G [^X ; 0]. Indeed, then, the distance between Dq = D 
and e* is at least Xi, and the distance between Dt and e* increases for t > (see Figure [6(b)] ). 

Let r be the smallest angle d can make with the plane containing Ti and refer to Fig- 



ure 6(c) Let At be the point of intersection between Dt and the plane containing Ti and Vt 
be the distance between At and the point on e* that realizes the distance between Dt and 
e*. The distance between Dt and e* satisfies d{Dt,e*) > wtsinF ^ c?(A(, e*) sinF. Hence, 
for proving that d{Dt,e*) ^ Xi for t ^ 0, it is sufficient to prove that d{At,e*) ^ for 
t ^0. We set a = to simplify the notation. 

We just proved that d{At, e*) ^ a implies d{Dt, e*) ^ Xi (for ^)- Conversely, we have 
that d{Dt,e*) < Xi imphes d{At,e*) < a. Similarly, for edge e, we get that d{Dt,e) < Xi 
implies d(^t,e) < a. By definition of Dt, we have that d{Dt,e) < Xi for ^ =^ 0, thus 
d{At,e) < a for i ^ 0. Furthermore, the angle between e and segment {At \ t E X jX"*"]} 
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(a) 



(b) 




Figure 7: For the proof of Lemma |4l 



is at least 7r/3 because this angle is at least the angle between their orthogonal projection on 
the plane defined by e and D that is the angle between e and D since all At lie in the plane 
spanned by Dt which projects on D] the lower bound of 7r/3 follows since the angle between 
e and D is at least 7r/3 by definition of Vl. Hence, the locus of points At, for t ^ 0, lies in 
a region, denoted T, shown in dark gray in Figure 7(a) For proving that d{At, e*) ^ a for 
t ^ 0, it is thus sufficient to prove that this region does not intersect the set, denoted T', of 
points at distance less than a from e 

and s be its 



Referring to Figures 7(b) (c) 



projection on the line supporting 



(shown in light gray in Figure 7(a) ). 
let p be the endpoint of e' the closest to e* 
e* . If the two regions T and T' intersect, there exists a 



point q in the intersection that is at distance less than or equal to -7=0 fromp and at distance 



V3' 

less than or equal to a from e*; thus d{p, s) ^ d{p, e*) ^ c?(p, q) + d{q, e*) ^ (1 



2 

V3- 



a. 



On 



the other hand, d{p, s) is one fourth of one of the heights of the triangle Tj and thus is at 
Hence, if the two regions intersect, then ^ ^1 



least 



2 

^/3 



Xi 

sinr ■ 



We postpone to 

Lemma [7] the proof that, with C = 31.3 (3, we have > ^1 + j i^^j which implies that 
the two regions T and T' are disjoint. This concludes the proof that no line Dt intersects 
the boundary of fn{Te)- 

We now prove that intersects fn{Te)- Consider a projection, p{), along the direction 
d onto a plane orthogonal to d. We proved that, for any of the two triangles Ti, T is at 
distance at least Xi from each edge e* 7^ e of . It follows that T lies in triangle Ti and thus 
that Dt intersects Ti for all t ^ 0. Therefore, D^ intersects Ti and is at distance at least Xi 
from each edge e* ^ e of Ti, for i = 1, 2. Furthermore, D~ is at distance x~ = miii(xi7 X2) 
from e, by definition. We now consider the triangle Ti for which Xi = X~ ■ It follows that D~ 
is at distance at least Xi from all three edges of Ti. Thus D~ projects to a point P — p^D") 
inside triangle p{Ti), at distance at least Xi from the three edges of p{Ti) (see Figure[8]). 
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Figure 8: Projection of dTi, fn{dTi) and D . 

Roughly speaking, by Hypothesis 2, the curve fn{dTi) is at distance less than Xi from 
dTi (the boundary of T^) tlius its projection p{fn{dTi)) is at distance less than Xj from the 
edges of p(Ti). It is thus intuitively clear that p{D~) intersects p{fn{Ti)), and thus that D~ 
intersects fn{Ti) (and thus fn{Tg)). 

More formally, consider the application Qn from the triangle p(Ti) to the plane containing 
it such that, for any point x in Tj, the point p{x) is sent to the point gn{p{x)) = p{fn{x)). 
We first prove that the curves p{dTi) and gn{p{dTi)) are homotopic in \ P. Consider the 
continuous application 

F : dTiX [0, 1] — > 

(x,A) Xp{x) + {l-X)gnip{x)) = Xp{x) + il-X)p{fn{x)). 

F is an homotopy between the curves p{dTi) and gn{p{dTi)) in M^. We prove that the image 
of F does not contain P, which yields the result. The triangle inequality gives 

d{P,F{x,X)) > d{P,p{x)) - d{F{x,X),p{x)). 

We have already proved that point P is at distance at least Xi from p{x) for all points x in 
dTi. On the other hand, the distance between p{x) and p{fn{x)) is larger than or equal to 
the distance between p{x) and their barycenter F{x, A), for any A G [0, 1]. Hence 

d(P, F(x, X))^x^~ d{p{x),p{,Ux))). 

Finally, since d{p{x),p{fn{x))) < Xi for all x G Tj, by Hypothesis 2, we have that, for all 
(a;, A), d{P, F{x, X)) > 0. Hence, the image of F does not contain point P and thus the 
curves p{dTi) and gn{p{dTi)) are homotopic in \ P. 

Now, we can contract p{dTi) to a point while remaining in p{Ti). Composing this with 
gn gives a contraction of gn{p{dTi)) in gn{p{Ti)). On the other hand, there is no contraction 
oi p{dTi) in \ P (since P is in p{Ti)), thus there is no contraction of its homotopic curve 
9n(p{dTi)) in \P. Hence, there exists a curve that is contractible in gn(p{Ti)) but not in 
IR2\P. It follows that gn{p{T^)) is not included in m?\P. Hence P is in gn{p{Ti)) = p{fn{Ti)). 
Therefore, D~ intersects fn{Ti) and thus fn{Te), which concludes the proof. 
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Figure 9: For the proof of Lemma [S] 



□ 

We finally prove the two following simple technical lemmas which complete the proof of 
Theorem [T] Recall that F is the smallest angle a direction d ^ can make with the plane 
containing Ti. 

Lemma 6. sin F = sin 77^= . 



Proof. In the following, we identify the sphere of directions with a sphere § embedded in 
R'^; let O denote its center. We assume that the embedding preserves directions (i.e., for 
any direction d, the corresponding point M on S is such that d and OM have the same 
direction) . 

Refer to Figure (5] Let d be a direction in 17 and M be its corresponding point on 
§. Consider one of the Ti and let 7 be the plane containing O and parallel to the plane 
containing Tj. Let H be the orthogonal projection of M onto plane T. Let E and E' be 
the two points on § that correspond to the two (opposite) directions of segment e. Let K 
be the orthogonal projection of M (and H) onto the line EE' . Finally, let be the angle 
ZMKH, (j) be the angle ZMOK, and 7 be the angle ZMOH. 

It follows from these definitions that 

HM HM KM . , . , 
sm 7 = — - — = — = sm sm m. 

' OM KM OM ^ 

Now, the angle 7 is also the angle between direction d and the plane that contains Ti. 
Thus F = inf 7, by definition of F. The angle 6 is the angle between the plane containing Ti 

and the plane containing e and d. It thus follows from the definition of that inf 9 = TrT= 



(see Figure 5(a) I. The angle (j) is the angle between d and the line containing e. It thus also 
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follows from the definition of fl that Jnf = f ■ In addition, since 7, 6 and (j> are in [0, 
we have 

C TT 

sinF = inf sin 7, inf sin 6* — sin ■— = and inf sin0 = sin — . 
rfeo rfGO 2yn den 3 

Furthermore, the constraints on 9 and (fi in the definition of O are independent. Thus, the 
minima of 9 and (t!) can be attained for the same direction d in Q. It follows that 



inf (sin 9 sin ch) ~ inf sin 9 . inf sin ( 

den desi den 



We can thus conclude that 



sin F = inf sin 7 = inf sin sin (p — inf sin inf sin (j) — —— sin ; 



den den den den 2 2^/n' 

Lemma 7. !f > (l + with C = 31.3/3. 

Proof. By LemmalHl replacing and F by their values in the inequality ^ > ^1 - 
gives 



□ 



V3 



^> (1 + 4 



§hi_ 

s/n 



4 V %/3; ^^sinf^ 

or equivalently 



15- sm , ^^-1= 



Notice first that for large enough values of rt, using the approximation sin a; « cc in the 
neighborhood of zero, we derive the sufficient condition 

Now, since we want our result for all rt, the computation is more complicated. Recall 
first that for any strictly concave function /, such that /(O) — 0, f(x) > "^^^ x for any 
X e (0, xo). It follows that sin a; > ^x for any x G (0, |). Since we assumed that 9e — -^ > 
and thus that < < % < -I , we get 

C \ 2 C 
sm I > 



2%/n ) IT 2yfn 
To guarantee inequality ([T|), it is thus sufficient to have 

4/3 1 + ^ ^k^f^- 



V3y 2 TT 2^' 
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or equivalently 




which concludes the proof. Note that we can set 



C ^ 31.3/3. 



(2) 



in the definition of il (in the proof of Lemma [3]) since Inequality llj is the only constraint 



3.3 Generalizations 

We prove here Theorem [2] We first show that Theorem [1] generalizes to the case where the 
viewpoint is chosen randomly at finite distance. We then show that considering surfaces 
with boundary does not change the asymptotic expected complexity of the silhouette. 

Point of view at finite distance. We have thus far restricted ourselves to the case where 
the viewpoint is chosen uniformly at random at infinity. However, our result applies to any 
distribution of viewpoints such that the probability for an edge e to be on the transparent 
silhouette is 0{9e), where 9e is the exterior dihedral angle associated to e; indeed, the 



when the point of view is chosen uniformly at random in a ball. This is also the case if S 
delimits a bounded region and the viewpoint is chosen uniformly at random in _B \ ^, for 
a ball B. 

Surfaces with boundary. Let S' be a 2-manifold with boundary H. We consider that the 
boundary is always on the transparent silhouette and so the definition of the transparent 
silhouette of a 2-manifold S with boundary is exactly that of a 2-manifold without boundary 
plus the boundary "B. 

The surface S is approximated by a triangulated mesh P„ that satisfies Hypotheses 1 
and 2, as in the case without boundary, except that now the mesh may not be a polyhedron 
(some edges may have only one adjacent face rather than two). 

To give an upper bound on the number of edges on the silhouette of the mesh, we consider 
the boundary edges and the other (non-boundary) edges separately. For the non-boundary 
edges, the same reasoning as before still holds. For the boundary edges, it is easy to see 
that the length (in 3D) of the boundary of P„ cannot be much larger than the length of "B. 
Indeed, the two are homeomorphic, and the hypotheses imply that the image of an edge e, 
of length le, is a curve whose extremities lie at distance at least 1^ — 2^- = apart. 

^Note that, in Lemma|3l silh(/n(Te)) always refers to an expected length for a viewpoint chosen randomly 
at infinity. 



on C. 



□ 




INRIA 



An Upper Bound on the Average Size of Silhouettes 



17 



This means that the length of 23 is at least r2(-^) times the number of boundary edges of 

Pn- Hence, the number of boundary edges of P„ is at most 0{y/n) times the length of 23. 
So, if the length of 23 is bounded, the expected complexity of the silhouette of P„ is 0{y/n). 

4 Conclusion 

This paper gives an idea of why, and when, the usual claim that the silhouette of a trian- 
gulated mesh has size 0{^/n) is valid. In particular, we have given a set of conditions such 
that any triangulated mesh approximating a surface in a way that satisfies those conditions 
has a silhouette of expected size 0(-\/n). Roughly speaking, the mesh should have no short 
edges, its faces should be fat, and the distance between it and the surface it approximates 
should never be too large. The surface itself is not necessarily everywhere differentiable and 
may have boundaries. 

A natural question to ask is whether meshes satisfying those conditions exist. In fact, for 
smooth surfaces, the meshes produced by Boissonnat and Oudot are one such example. 
The critical property of the meshes they compute is that the ratio between the size of the 
largest and the smallest triangles remains bounded, although meshes are non-uniform with 
small triangles in areas of large curvature. However, in order to satisfy our conditions, non- 
smooth surfaces with curved sharp edges (such as a flying saucer with a sharp equatorial arc) 
would have to be approximated by small triangles over the whole surface. Such meshes would 
have silhouettes of expected size 0{\/n) but then n would be much larger than necessary; it 
would be reasonable to replace the large number of triangles used to mesh large flat portions 
of the surface with a smaller number of large triangles, which would give a silhouette of size 
closer to linear. This explains why the observed expected size of silhouettes, as shown in 
0, is larger than 0{^/n). The fact that non-uniform meshes approximating such surfaces 
appear, in computer graphics, to have silhouettes of expected size much smaller than n is 
thus likely due to additional properties of the surfaces or the meshes. 
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